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^■. ABSTRACT 
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^ ' Gamma-ray burst (GRB) afterglow observations have allowed us to signifi- 

p ■ cantly constrain the engines producing these energetic explosions. Te redshift 
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-<^ ■ and position information provided by these afterglows have already allowed us 

O ' to limit the progenitors of GRBs to only a few models. The afterglows may 

also provide another observation that can place further constraints on the GRB 
progenitor: measurements telling us about the environments surrounding GRBs. 
pi ' Current analyses of GRB afterglows suggest that roughly half of long-duration 

^2 ' gamma-ray bursts occur in surroundings with density profiles that are uniform. 

t::^- ■ We study the constraints placed by this observation on both the classic "col- 
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lapsar" massive star progenitor and its relative, the "helium-merger" progenitor. 



O ■ We study several aspects of wind mass-loss and find that our modifications to the 

^ ■ standard Wolf-Rayet mass-loss paradigm are not sufficient to produce constant 

I ■ density profiles. Although this does not rule out the standard "coUapsar" progen- 

O 

c/3 . models and find that such progenitors can fit this particular constraint well. We 

show how detailed observations can not only determine the correct progenitor for 



itor, it does suggest a deficiency with this model. We then focus on the He-merger 



GRBs, but also allow us to study binary evolution physics. 

Subject headings: binaries: close — black hole physics — gamma rays: bursts- 
stars: neutron 



1. Introduction 

Since the concurrent and cospatial observation of supernova 1998bw with gamma-ray 
burst 980425 (Galama et al. 1998), evidence that long-duration gamma-ray bursts (GRBs) 
arise from the explosion of massive stars has continued to grow. This supporting evidence 
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includes the association of long-duration GRBs with star-forming galaxies (Fruchter et al. 
2002; Vreeswijk et al. 2001) and the detection of 1998bw-like "bumps" found in GRB after- 
glows (Bloom 2003). The connection between long-duration GRBs and massive stars was 
sealed by the paired observation of SN 2003dh and GRB 030329 (Hjorth et al. 2003; Stanek 
et al. 2003). This association tells us that the progenitors of long-duration GRBs are al- 
most certainly massive stars. Both theory and observation further constrain the nature of 
these massive star GRB progenitors. A successful GRB engine must produce beamed jets 
of relativistic matter with Lorentz factors in excess of 100 (Frail et al. 2001). To reach 
and maintain the high Lorentz factors observed in GRB ejecta, the jet must not sweep up 
much mass and the jet engine must be active for the entire time during which the jet pushes 
through the star, leading theorists to argue that the massive star must lose its hydrogen 
mantle before exploding (Woosley 1993; Zhang, Woosley, & McFadyen 2003). Such a star 
would have been a Wolf-Rayet star before collapse and its resultant supernova explosion 
should have matched type Ib/c supernovae. Spectral observations of SN 2003dh have con- 
firmed that the progenitor was a Wolf-Rayet progenitor (Hjorth et al. 2003; Stanek et al. 
2003). 

Wolf-Rayet stars are known to have strong winds with mass outflow rates beyond 
lO~^M0y~^ and wind velocities above 1000km s~^. Stellar winds sweep out the environ- 
ment surrounding the star and, depending upon the density of the ambient medium, these 
wind-blown bubbles can extend beyond 10 pc (Dwarkadas 2005). The density profile for 
a free-streaming wind is simply determined by mass conservation: M„ind = 47rr^rp(r) -^ 
pij-) = M„ind/47rr^t>wind oc r~^ where Mwind is the mass loss rate of the wind (assumed to 
be constant) and p, f wind are the density and velocity (or r) of the wind respectively. The 
density profile takes on a r^^ radial dependence when we assume a constant wind velocity. 

Very quickly, GRB fireball theorists (see Chevalier et al. 2004 for a review) realized 
that this density profile, with its characteristic p oc r^^ profile, would produce a different 
afterglow light curve than, for example, a constant density profile expected from alternate 
GRB models (e.g. the merger of a double neutron star system). With all of the evidence 
pushing toward massive star progenitors for long-duration GRBs, we would expect to find 
that GRB afterglows were better fit by p ex r~^ density profiles than by constant density 
profiles, but roughly half of all systems seem to be better fit by constant density profiles, with 
typical densities between 0.1 and 100cm~^(Meszaros et al. 1998; Dai & Lu 1998; Chevalier 
& Li 1999; Panaitescu & Kumar 2002; Chevalier et al. 2004). The disagreement between our 
picture of massive star progenitors for GRBs and observations of the GRB afterglow remains 
an outstanding problem in our understanding of long-duration GRBs. In this paper, we will 
study predictions of the environments surrounding long-duration GRB progenitors in an 
effort to resolve the density-profile problem. 
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A number of different approaches might resolve this particular discrepancy. On one 
extreme, one can ignore the observations, hoping that future data will show a different trend 
and negate the current results^. On the other extreme, one can rule out the current paradigm 
for long-duration gamma-ray bursts and search for a different, non-stellar progenitor. In 
between a number of options exist and we will pursue two of these: 1) our free-streaming 
p oc r~^ profile estimate for the density in the wind-blown region is an oversimplification; 
or 2) the progenitor does not reside in its wind-blown bubble when the GRB is produced. 
We further restrict our study to progenitors of the collapsar or collapsar-like GRB engine 
(Woosley 1993; Fryer & Woosley 1998; MacFadyen & Woosley 1999). This engine, which 
invokes the collapse of a massive star down to a black hole and the formation of an accretion 
disk around that black hole, only works for a small subset of massive star progenitors. We 
believe that by studying this small subset of progenitors, we may gain some insight into the 
environments surrounding GRB engines. 

The first option, modifying the density profile in the wind-blown bubble, can be ac- 
complished in a number of different ways and has already been studied at some level. Some 
observations can be made to fit a wind profile if the density is decreased (Wu et al. 2003; 
Dai & Wu 2003). Chevalier et al. (2004) have also, in a generic sense, studied the true 
termination radius of the Wolf-Rayet wind. If the wind terminates at a small radius, the 
observations are probing the constant density region beyond the termination shock, and not 
the p oc r~^ free-streaming wind region. In §2 of this paper, we study the wind profiles of 
massive stars designed to produce coUapsars ("massive star" progenitor to the collapsar en- 
gine), using a combination of stellar evolution models and 1-dimensional and 3-dimensional 
wind models. The wind mass loss is time-dependent, and we study the effects of this time 
dependence. Most working collapsar models require binary interactions, and we also study 
these binary effects. 

Alternatively, there is a collapsar progenitor scenario where, before the GRB explosion, 
the progenitor moves out of its wind nebula: this is the helium-merger progenitor (Fryer & 
Woosley 1998; Zhang & Fryer 2001). In this progenitor, the more massive star in a binary 
collapses to form a neutron star or black hole, providing a kick to its binary system. Not 
until its companion evolves off the main sequence and causes the stars to merge is a burst 
produced. Although Fryer & Woosley (1998) already predicted that such a progenitor would 



^A number of assumptions go into the afterglow models needed to estimate the density profiles and it 
could well be that these assumptions are incorrect. Although roughly half of all explosions are better fit 
by a wind density profile {p oc r^^), constant density profiles can not be ruled out for any gamma-ray 
burst. But for Only 1 or 2 cases can a r^^ density profile be ruled out emphatically (Panaitescu - private 
communication). This observation is far from certain. 
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more likely yield a surrounding environment akin to a constant density interstellar medium 
(ISM) environment, no work to date has studied this claim in detail. In §3, we show results 
of Monte Carlo calculations to study the positions of these bursts relative to the formation 
site of their binary progenitor. 

We ultimately find that the seeming discrepancy between simple wind models and GRB 
afterglow observations can play an important role in homing in on the true progenitor of 
these long-duration bursts. Although we can not rule out the massive star progenitor GRB 
scenario, it is becoming increasingly difficult to match such a progenitor to observations. 
The helium merger scenario is best able to fit current afterglow lightcurves, but it may have 
problems matching other observational constraints. We conclude by discussing the issues 
remaining for the different long-duration GRB scenarios. 



2. Stellar Winds from Collapsar Progenitors 

Weaver et al. (1977) first discussed the structure of wind-driven bubbles produced by 
massive stars. Assuming that the wind profile varies with at most a linear dependence on 
time and is spherically symmetric, the structure of these wind-blown bubbles is well known 
(Koo & McKee 1992a, 1992b). However, as we shall see in this section, neither of these 
assumptions is strictly true. We study the bubble structure produced by the rapidly varying 
winds from Wolf-Rayet stars, the likely progenitor star for collapsar GRBs. 



2.1. Understanding Wind Profiles 

Before we study realistic wind profiles based on stellar mass-loss histories, we first build 
our intuition of bubbles from winds with constant mass-loss rates blowing into environments 
with different densities. Such wind blown bubbles have been studied in detail analytically 
(Koo & McKee 1992a, 1992b; see Dwarkadas 2005 for a review) and the wind environment 
can be described by 4 separate regions: (1) the freely- flowing fast wind, (2) the shocked 
fast wind, (3) the shocked interstellar medium and (4) the unshocked interstellar medium. 
For the purpose of comparing to the afterglow observations, the primary constraints come 
from the first region, where the free-fiowing winds develop a p oc r~^ density profile. To 
overcome the observational constraints, we must either lower the density of this region or 
limit its size. The free-fiowing region is capped by a termination shock (rtermination), and 
we will focus on determining what conditions can drive down this termination shock to low 
radii: rtermination < 0.1 pc. At such small size scales, the GRB jet can pass through this 
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region without seriously impacting the afterglow emission. 

To understand the dependence of the termination shock radius on the uncertainties in 
wind parameters, we have modified the 1-D Lagrangian code developed by Fryer et al. (1999) 
to include the addition of mass from winds. Wind material is added in zones of constant 
mass. For constant mass loss from winds, this means that zones are added at constant 
time intervals. The inner boundary of the simulation moves out with the wind velocity. 
When a particle is added, its outer radius is set to the inner boundary radius and the inner 
boundary is set to the initial inner boundary radius. By varying the time intervals according 
to individual wind velocities, this technique can be used to simulate winds with varying 
mass loss rates and varying velocities. We model a region from ~0.01 to 20-200 pc. For our 
1-dimensional, constant mass-loss studies, the interstellar medium is modeled with between 
~50,000 and 100,000 zones. In 250,000 y, roughly 25,000 zones are added to the simulation. 
Table 1 summarizes all of the parameter study models. 

Figure 1 (top) shows the density profiles for four different wind bubble simulations 
differing only in the density of the ambient medium (roughly 10, 100, 1000, and 10,000 
cm^^). We have assumed a constant wind velocity (1000 km s^^) and constant mass loss 
rate (10~^Mq|/~^) which roughly corresponds to what we would expect from a Wolf-Rayet 
wind. A lot can be learned from these simulations, but we will focus on the radius of 
the termination shock (rtermination) • 250,000 y of this strong wind produces a termination 
shock that is at a radius of 1 pc if the wind is blowing into an interstellar medium with a 
density of 100 cm~^. Fitting to these simulations, we find that the radius of this termination 
shock is ~ l(nisM/100cm~^)~^/^pc. It would require very high densities indeed to place the 
termination shock below 0.1 pc. We include a final simulation (dotted line) that includes a 
cooling routine outlined in Fryer et al. (2006) based on the cooling rates of Sutherland & 
Dopita (1993). The primary effect of cooling is to alter the density profile near the shock 
with the interstellar medium. Another effect is the fact that the UV photons from the star 
will drive an ionizing front ahead of the wind. This has been recently studied by van Marie 
et al. (2004). Although much more must be done to understand its effects, van Marie et al. 
(2004) found that the primary change from the ionization front was in the outer radius of 
the wind shock. It does not significantly alter the radius of the free-streaming region and we 
will not consider it further in this paper. 

Figure 1 (bottom) shows the dependence of the termination shock radius on the other 
free parameters in a constant wind: the mass-loss rate (solid lines) and the wind velocity 
(dotted lines). Here we have assumed a low density interstellar medium (10cm~^). The 
dark solid line assumes a mass-loss rate of 1O~^M0 y^^ whereas the light solid line assumes a 
mass-loss rate of lO'^M© y~^. The termination shock radius is also proportional to the mass 
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loss to the 1/3 power. The dark and hght dotted hnes assume wind velocities of SOOkrns"^ 
and 2000 km s^^ respectively. Although this changes the position of the outer shock of the 
wind-blown bubble, it does not significantly affect the radius of the termination shock. 

To restrict the termination shock to a radius below 0.1 pc, either the ambient density 
surrounding the GRB progenitor must be high {riisu > 10, OOOcm^^) or the Wolf-Rayet wind 
must have a mass- loss rate much lower than the canonical IO^^MqI/^^. Such constraints 
may not be that extreme. Densities in excess of 10, OOOcm^^ are common in molecular 
clouds, especially near massive star birthsites. For example, the shock producing the young 
supernova remnant Gas A is impacting a 650 M0 molecular cloud with densities of lO^cm"^. 
Unfortunately, such a dense medium is ruled out already by afterglow observations (Berger 
et al. 2003; Chevalier & Li 2000; Panaitescu & Kumar 2002; Kobayashi & Zhang 2003)^. 

However, the Wolf-Rayet mass-loss rate could be much lower than the canonical value 
for the specific stars that make GRB progenitors. Such a low mass-loss rate might result 
from a helium star formed in a binary interaction that causes a 15-20 Mq star to become a 
Wolf-Rayet star. Such a star can have a Wolf-Rayet mass loss rate as low as 10^"^ MQy^^. 
This is not entirely satisfactory, as we believe such a star is not likely to collapse to a black 
hole. Alternatively, because the Wolf-Rayet mass-loss rate depends upon the metallicity 
(Mwoifflayct oc Z^'^ — Z^'^), GRB progenitors with observed surroundings matching constant 
density profiles could be limited to low metallicity stars (Yoon & Langer 2005). A low 
metallicity star turned Wolf-Rayet via a binary common envelope phase could have a very low 
Wolf-Rayet mass loss rate. For example, if Mwoif-Rayct c<: (Z)°-^ (where Z is the metallicity), 
then a 25 M0 star with a metallicity 1/lOOth that of solar could have mass loss rates below 
IO^^Mq^/"^. Not only will the free-streaming region of this wind be small, but its total 
density will also be low (a better fit to the observations). With such stars, ambient densities 
of a few hundred cm^^ would be sufficient to restrict the termination shock to less than 
O.lpc. 

Before we move onto the complexities arising from realistic wind models, we mention 
a few other effects that may also limit the free-streaming wind profile in GRB progenitors. 
These stars may only form in very dense clusters and, if so, the interaction of the winds 
could limit the extent of the free-streaming profile. Rockefeller et al. (2005) studied the 
hydrodynamics of winds in the Arches and Quintuplet clusters. By analyzing their data, we 
find that the free-streaming wind region of even the strongest winds can be capped at 0.5 pc 



^Although bear in mind that if we are observing the GRB shock as it progresses through shocked wind 
and ISM material, the observed density wiU be much less than the density of the interstellar medium. Thus, 
it is possible that large interstellar medium densities can exit. 
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at the center of the cluster. 

In addition, these massive stars may have large velocities with respect to their circum- 
stellar environments. The ram pressure of the wind environment against this circumstellar 
environment limits the extent of the wind. This extent can easily be derived by setting the 
free-streaming wind pressure = O.Spwind'wiind = ^wind/(87r)r^j^^-i;„ind to the ram pressure 
of the motion through the circumstellar medium = O.SpisM'^star- Here Pwind,'ywind!''"wind and 
Mjvind are, respectively, the density, velocity and radial extent and mass-loss rate of the free- 
streaming wind; pisM is the density of the interstellar medium, and f star is the motion of the 
star with respect to this circumstellar medium. Setting these two ram pressures equal, we 
derive an equation for the radial extent of the wind: 

M„ind \ ( t^wind ^ ■ ( PiSM \ ~°-^ ( ^star \ ~^ ., ^ 

lO-SM^y-i; VlOOOkms-V VlOcm-3; UoOkms-iy P"' ^'^ 

The number of massive stars with velocities in excess of lOOkms^^ is very low (10 — 30%: 
Hoogerwerf et al. 2001) and unless something in the progenitor evolution requires these fast 
velocities, it is unlikely that these motions can explain the environments. 



2.2. Real Stellar Mass Loss 

With the intuition gained from our constant mass-loss winds, we are set to attack the 
more realistic problem of stellar winds from GRB progenitors. In this paper, we use the 
output from the Tycho stellar evolution code (Young & Arnett 2005) and we study three 
different possible GRB progenitors: a 40 Mq single star, a 23 M© binary system and a 16 Mq 
binary system. In the binary systems, we assume that a common envelope phase removes the 
hydrogen envelope of the primary when the primary star expands off the main- sequence. The 
mass loss rate and wind velocities for the GRB progenitors are calculated at each timestep 
during the evolution. Three prescriptions are used during steady mass loss phases. Mass 
loss for blue, hydrogen-rich stars {logT^ff > 3.875) uses Kudritzki et al. (1989). At lower 
effective temperatures appropriate for the red giant /supergiant phases, the prescription of 
B15cker (1995) is used. For Wolf-Rayet mass-loss rates, we use the Lamers & Nugis (2003) 
prescription. 

Catastrophic mass loss episodes are dealt with in a more schematic fashion. The 40 Mq 
star is subject to luminous blue variable eruptions. The conditions for the instability (large 
radiative accelerations at temperatures coinciding with "bumps" in the continuum opacity) 
are idenitifiable in the stellar model. When these conditions are satisfied a large amount of 
mass is removed from the model to mimic an eruption. Limitations on the numerical stability 



of the code prevent the mass from being removed as quickly as in a real luminous blue variable 
(LBV) eruption. Instead mass is lost at a high but steady rate dictated by code stability 
for a period of order 100 years with exact duration dictated by the amount of mass loss 
necessary. Such episodes can be identified in Figure 2. We choose a conservative approach, 
removing only enough mass to stabilize the stellar structure, though it is possible that the 
envelope could be removed to the depth of the instability generating region. Depending on 
the stellar structure, this could represent several solar masses of ejection. In the first, large 
eruption ~ IMq is removed. Subsequent eruptions are smaller, and the bulk of the mass 
loss occurs via steady winds. 

The 16 and 20 Mq models are intended to represent stars which lose their hydrogen 
envelopes during the first ascent Red Giant Branch (RGB) to a binary interaction. As with 
catastrophic eruptions, the physical processes cannot be captured realistically in a ID stellar 
evolution code, so we must use an approximation. A similar procedure is used, which results 
in an ejection of the hydrogen-rich envelope (~ 10 Mq , depending on stellar mass and mass 
loss history) in less than a thousand years. From an evolutionary standpoint, this is a good 
approximation, as the nuclear timescale of the star is much greater than a thousand years 
at this stage. This is not necessarily a good approximation to the mass loss, but it does 
provide the correct qualitative behavior: a short period of very high density, low velocity 
(~ Vescape) uiass loss which clcars a cavity around the star into which a fast Wolf-Rayet 
type wind expands. The star subsequently evolves as a WR towards a Type lb supernova. 
Though the binary ejection of the hydrogen envelope is handled only in a very schematic 
way it does represent, along with the subsequent evolution, a significant improvement over 
constant density or power law density profiles for the circumstellar medium. 

The Wolf-Rayet mass loss is characterized by a high mass loss rate, high velocities, and 
an episodic history. The abrupt changes in mass loss occur when the main driver for the mass 
loss switches between opacity bumps at different temperatures and hence different radii in 
the star (see Lamers & Nugis 2003 for a discussion of the Wolf-Rayet mass loss mechanism). 
Early transient spikes may represent eruptions driven by high radiative accelerations in the 
stellar atmosphere akin to LBV or Ofpe (Massey et al. 2000) eruptions. 



2.3. Wind-Blown Bubbles From Real Winds 

The winds from our Tycho models have a number of characteristics that will complicate 
the wind profile. Note that the ejection velocity is often low when the mass loss peaks (either 
from a LBV eruption or a common envelop mass ejection). The high- velocity wind ejecta 
that follows these peaks will bounce against the slow-moving, dense, peak ejecta, causing a 
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shock to flow back down through the wind. Also, a similar shock will be produced when the 
eruption hits the termination shock. With realistic winds, our free-streaming, "unshocked" 
wind region now has many shocks running through it. This will dramatically affect the 
entropy profile of this wind. But afterglow observations are not sensitive to the entropy, but 
rather the the density profile of the wind-blown bubble region. 

To study the density profile, we first return to our 1-dimensional simulations. By al- 
lowing non-constant time intervals for the addition of zones, we can incorporate the variable 
mass-loss rate and time-dependent wind velocity profiles from our Tycho models. We assume 
an ambient interstellar density of lOOcm"^. Table 1 lists the simulations of these realistic 
mass-loss profiles; we give rough ranges for the mass loss rates and velocities instead of ab- 
solute values. We model the wind-driven bubble only during the last few hundred thousand 
years of the star's life: 0.6, 0.5, 0.6 Myr for the 16,23,40 Mq stars respectively. Although 
this reduced time does not accurately model the outer extent of the winds, it does get the 
position of the free-streaming region fairly accurately^ . Most of the explosive phases of mass 
loss occurs during this time, and it is this ejecta that dominates the density profile of the 
free-streaming wind region. Figure 3 shows the density profiles of these models fiowing in 
an ambient density of 100 cm"'^. For these models, it is difficult to pinpoint the termination 
shock of the free-streaming wind and there is no strict p oc r~^ dependence to the density. 
But the density does decrease roughly as r~^ well beyond 1 pc. Recall that the termination 
shock from our 1Q~^ MqU'^^tii^m = 100 cm^^ calculation was at Ipc (even for the 40 M© 
star where the mass loss exceeds this standard value). The ambient medium would have to 
exceed lO^cm"^ to reduce this region below O.lpc. This density is probably ruled out by 
the timescale of the radio emission (see Section 2.1). 

We find that episodic spikes in the mass loss produce a termination shock which causes a 
significant departure from an p oc r^^ density profile. This does not help the situation much 
because the region of departure from a free-streaming profile propagates out from the star at 
roughly the velocity of the slow moving ejecta (~ a few 100km s~^). Thus the density profile 
is disturbed within the O.lpc limit for 1000 years or less, roughly 1% of the time between 
eruptions. Although we cannot rule out such epsiodes within 1000 years of collapse unless 
the stellar surface is hotter than ~ 5 x 10^ K (hotter than the continuum opacity bumps 
which drive the high mass loss rates), the odds of an ejection event occurring within 1000 
years of collapse are too small to explain the distribution of afterglows unless the frequency 
of ejection events increases as the star approaches collapse. We do not observe this behavior. 



■^Vikram Dwarkadas is currently modeling the full wind profiles of our binary star in an effort to under- 
stand the supernova remnant Cassiopeia A. The extent of the free-streaming region from his shock agrees 
with ours at the 20% level. 
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but the uncertainties in the physics of mass loss do not allow us to eliminate the possibility. 

Another wrinkle in this model is the fact that these winds are far from spherically 
symmetric. As argued by Fryer et al. (1998), the idea that GRB progenitors most likely 
arise from binary systems is gaining increasing acceptance (Petrovic et al. 2005, Fryer & 
Heger 2005). If the wind is not already asymmetric due to rotation (Maeder & Meynet 
2000), the binary interaction will ensure that it is. If the mass ejection arises from common 
envelope evolution, then this mass outflow at least is highly asymmetric. LBV eruptions are 
also observed to be strongly bipolar. 

We have modified the three-dimensional SNSPH code (Fryer, Rockefeller & Warren 
2005) to model these winds and then study the effects of asymmetries on the density struc- 
ture. We add particles in shells centered on the location of the wind source, after a time 
At ~ h/vwind has passed since the last addition of particles, where h is the SPH particle 
smoothing length and Vwind is the current value of the wind velocity. The mass of each par- 
ticle is set according to m = M^ind^t/nsheii, where M^ind is the current value of the wind 
mass loss rate and ngheii is the number of particles per shell. We remove particles from the 
simulation when they pass beyond 3 pc from the wind source. 

Figure 4 shows the density profile of one such SPH simulation. A 23 Mq binary system, 
located at the center of the simulated volume, ejects material during a common envelope 
phase at a rate of ~ O.IMq yr^^. We assume that this material is not ejected isotropically; 
instead, we place it in a 45° wedge oriented perpendicular to the z axis. The subsequent 
wind from the binary system has a much higher velocity but represents a much lower rate 
of mass loss from the system, so the wind collides with the inner edge of the dense wedge 
and flows out of the system primarily along the z axis. In this calculation, each added shell 
of wind material contains 400 SPH particles, each shell is added at a radius of 0.01 pc from 
the wind source, and the inner edge of the dense wedge is at r = 0.3 pc. The total number 
of particles in the simulated volume is ~ 350, 000 but fluctuates slightly as the wind velocity 
changes. 

Despite the modifications we have introduced to the wind profiles, we have not produced 
any conditions that cause the density profile of the wind-blown bubble to differ significantly 
from the p oc r^^ dependence. With these results in mind, our proposed solution for a 
the massive star progenitor for GRBs is a low-metallicity, 20-25 Mq star turned into a Wolf- 
Rayet star through a binary common envelope phase. A 23 Mq progenitor with a metallicity 
1/lOOth that of solar would have a weak Wolf-Rayet wind, but also be massive enough to 
collapse to form a black hole. A GRB in a high metallicity host would rule out such a 
scenario, forcing us to progenitor models like the He-merger scenario. 
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3. He-Merger GRB Progenitors 

Before we study the environments surrounding He-merger GRB progenitors, we first 
review the He-merger scenario itself. The He-merger progenitor for GRBs invokes a binary 
system where the most massive star collapses to either a neutron star or a black hole. If this 
compact object merges with its companion, it will spiral into the center of the companion 
while spinning up the companion. In this merger, the binary has satisfied the three main 
requirements behind the collapsar model: 1) the compact object accretes enough to collapse 
to a black hole, 2) the star surrounding the black hole has enough angular momentum 
(perhaps too much as we will discuss below) to form a disk around that black hole, and 3) 
the hydrogen envelope is ejected in a common envelope excretion disk. 

There are two main binary evolution scenarios that dominate the production of He- 
Merger progenitors of GRBs: common envelope mergers and kick-induced mergers of nearly 
equal mass stars. Figure 5 shows the basic evolutionary phases in the common envelope 
scenario. In this scenario, we generally begin with an initial binary in a fairly close orbit (< 
1 A.U.). The more massive star ("primary") evolves off the main sequence and, for most 
progenitors, envelops its companion as it turns into a giant (or supergiant) star. This causes 
the orbital separation of the binary to tighten in a common envelope phase. The primary 
then collapses, forming either a neutron star or black hole. A fraction of these binary systems 
will remain bound after the supernova explosion. When the companion star evolves off the 
main sequence, it can envelop the compact remnant produced by the collapse of the primary. 
A fraction of these binaries will merge during this second common envelope phase, producing 
a helium- merger GRB. 

The second main scenario we loosely term the "Brown" He-merger scenario because it 
uses similar initial conditions (a binary system with nearly equal massed stars) to a double 
neutron star formation scenario proposed by Bethe & Brown (1998). The idea of Bethe & 
Brown (1998) was to consider binary systems where the two stellar components have masses 
that are so close that the less massive star evolves off the main sequence before the more 
massive star collapses to form a compact remnant (this corresponds to masses of the two 
stars that are within ~10% of each other). In their model, the two helium stars produced 
during the mass transfer phase are pushed into a tight binary that remains bound through 
the two subsequent supernova explosions, forming a system with two compact objects. In 
our scenario, we follow the Bethe & Brown (1998) evolution up to the point that we form 
two helium stars in a tight binary (Fig. 6). When the first helium core collapses, its kick 
can drive it into its companion, producing a helium-merger GRB. 

This "Brown" He-merger scenario uses supernova kicks to cause the compact object to 
merge with its stellar companion. This begs the question of whether a kick imparted to the 
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compact object formed during the collapse of the primary can drive the neutron star or black 
hole into its main sequence companion (Fig. 7). Because main sequence stars are not quite 
as compact as helium stars, we expect lower accretion rates in these models, possibly making 
this scenario less viable as a GRB progenitor. In addition. Fryer et al. (1999) pointed out 
that this main-sequence star mechanism is rare. We will explicitly discuss the rates here, 
finding that this mechanism accounts for less than ~5% of the total helium-merger GRB 
population for most simulations. 

These three different progenitor scenarios produce very different environments through 
which the GRB jet will pass (Table 2). Both mechanisms 2 and 3 produce mergers very 
soon after (generally within a few years of) the supernova explosion. GRB explosions for 
these two scenarios will first pass through a young supernova remnant and then through the 
wind from the Wolf-Rayet progenitor. Scenario 1 does not produce a GRB until long after 
(~ 1 — 10 Myr) the collapse of the primary star. With the kick imparted to the system during 
the collapse of the binary, most bursts from this progenitor occur beyond the edge of the 
supernova remnant and the Wolf-Rayet wind termination shock. Hence, the environment 
surrounding these bursts will not be the wind-swept structure characteristics we expect 
from Wolf-Rayet stars. The immediate environment surrounding these bursts will either be 
characterized by a constant density medium (if the companion is low-mass: < 15 — 20 M©) 
or by the giant winds from the companion. In this section, we study the relative frequencies 
of occurrence of these different burst scenarios and their velocity, merger time, and position 
distributions including a brief study of the uncertainties of all these values. 



3.1. Possible Constraints on the Helium merger 

It is likely that only a subclass of He-mergers will produce GRBs. Currently we do not 
understand the progenitors or the collapsar mechanism well enough to definitively exclude 
some progenitors. Perhaps observations of GRBs can provide some insight into the progenitor 
subclass that truly makes GRBs. Here we discuss a variety of He-merger subclasses that 
may make up the true set of GRB progenitors. We will then study each of these subclasses 
separately when we present our population synthesis results. Table 3 summarizes these 
different subclasses. 

One of the sources of uncertainty in our population synthesis studies is that we don't 
know exactly what conditions actually produce GRBs. For example. Fryer & Woosley (1998) 
assumed that the primary could collapse into either a neutron star or a black hole and still 
be a candidate GRB progenitor. Their assumption, which has been borne out by current 
simulations (Zhang & Fryer 2001), is that even if the compact remnant were a neutron star. 
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it would accrete so much in the inspiral phase that it would become a black hole quickly, 
satisfying the "black hole requirement" producing coUapsars. But if the accretion rate is not 
so extreme (angular momentum, for example, could drive outflows), then the neutron star 
may not collapse to form a black hole and such a progenitor will not produce a collapsar. 
Alternatively, it may be that black holes accrete too much mass, making inefficient collapsars 
unable to drive strong GRB jets. We will look separately at systems whose primary collapses 
to a neutron star versus those systems whose primary collapses to form a black hole. This 
gives us two separate subclasses based on the initial status of the compact remnant: 1) a 
neutron star or 2) a black hole. 

Another constraint on the progenitor may be the amount of angular momentum. When 
Zhang & Fryer (2001) modeled their mergers, they focused on maximizing the angular mo- 
mentum surrounding their collapsar black hole. They assumed that the stars were, if any- 
thing, spinning in the same direction as the orbital angular momentum. They also assumed 
that very little orbital angular momentum was required to remove the hydrogen envelope 
and that all of the orbital angular momentum was deposited in the helium star. This pro- 
duced very high angular momentum disks in the merged black hole accretion disk system, 
leading some to argue that this progenitor does not produce a viable GRB (Di Matteo et 
al. 2002). However, recall that Zhang & Fryer (2001) were trying to maximize the angular 
momentum in the system. The angular momentum could be much lower if most of the or- 
bital angular momentum is lost to the hydrogen envelope or if the star is spinning in the 
opposite direction with respect to the binary orbit. Mergers produced with kicks may also 
lead to lower angular momenta. In general, we assume that a star's spin direction in a binary 
coincides with the direction of the orbit. This is almost certainly true prior to the supernova 
explosion. But the kick imparted to the compact object during supernova explosion can (and 
often does for those binaries that remain bound) produce binaries where the spin angular 
momentum direction is 180 degrees off from the orbital angular momentum of the binary. 
This subclass is defined in our calculations by restricting our sample to those companion 
stars whose helium star moment of inertia (from the helium stars of Fryer & Heger 2005, 
-^He-Star ~ O.lMHe-star-Rnc-star) ^^ greater than the moment of inertia of the merging compact 
object when it reaches the helium star radius (/compactBinary = Mcompactobjcct-Rnc-star)- 

Lastly, if the companion star is too massive, it too will blow a wind, leading to a 
surrounding wind environment not too different from the Wolf-Rayet winds produced by the 
primary star. We define this subclass of low-mass losing systems by restricting our sample 
to systems with companions whose total mass-loss is less than 1 Mq. This essentially limits 
us to binaries where the initial mass of the companion is less than 15 M©. 

With these four subclasses in mind, we can now analyze our population synthesis studies. 
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3.2. Population Synthesis studies 

Our population synthesis calculations are based on the code described in Fryer et al. 
(1998,1999). The main alteration to this code has been to update the mass of the collapsed 
core using the results of Fryer & Kalogera (2001). Fryer et al. (1999) set the remnant mass 
for stars collapsing to neutron stars to 1.4 Mq and the remnant mass for stars collapsing to 
black holes to one third the mass of the star at collapse. In this paper, we have fit a function 
to the results of Fryer & Kalogera (2001), allowing a more continuous spread of neutron star 
and black hole masses. We refer the reader to Fryer et al. (1998,1999) for any other details 
about the population synthesis technique or the specific code we use. We have run over a 
dozen models using a wide range of choices for the uncertain binary population synthesis 
parameters (Table 4). 

Table 4 shows a list of all the population synthesis models studied in this paper followed 
by the values of the population synthesis parameters used by each model. For a more 
detailed description of each of these parameters, see Fryer et al. (1999). We focus on just 
a few parameters: the power law of the initial mass function (aiMp), the distribution of 
velocities imparted to compact objects at collapse (kick), mass transfer parameters denoting 
the fraction of overfiowing matter in a binary accreted by the companion (/5mt) and the 
specific angular momentum of the overfiowing matter ejected, and hence lost, from the 
binary (c^mt), and the efficiency with which the inspiralling star in a common envelope 
phase ejects the envelope (q;ce)- We sample both our primary and secondary stars from the 
same initial mass function, constraining the secondary to be less massive than the primary. 
After some argument about the distribution of kicks arising from stellar collapse, it is now 
generally accepted that a bimodal kick distribution like that proposed by Fryer, Burrows, 
& Benz (1998) is roughly correct (Arzoumanian et al. 2002). We term the kick distribution 
represented by two Maxwellians of 50 and 500 km s^^ the FBB kick distribution. We have 
also used a set of single Maxwellian distributions covering a range of mean values to better 
understand the of the dependence on the kick magnitude. Finally, we have run a few tests 
of the dependencies on the stellar mass loss rates and radii. We model 50 million binaries in 
each simulation. 

Fryer et al. (1999) presented very few results outlining the characteristics of the Helium 
Merger model, focusing mostly on the formation rate. Here we study the distribution of 
merger velocities, merger times and distances traveled before GRB outburst for the various 
subclasses of He-merger progenitor. The explosion and the supernova kick impart a velocity 
to the binary (if it remains bound). Figure 8a shows the distribution of the post-collapse 
velocities of He- Mergers for our "Standard" model (Table 4). The dark solid line shows 
the number of binaries per Ikms^^ bin for all He-merger progenitors. Although for most 
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subclasses, the bulk of the systems have velocities below 25kms~^ (where the total number 
peaks), a second peak exists above lOOkms"^. The two peaks arise because of the double 
peaked nature of the supernova kick distribution. The peak at low velocities is not present for 
two subclasses of He-merger models: those with initially black hole compact remnants and 
those with high moments of inertia. Note that the peak in low- velocity systems is dominated 
by the NS Collapse and Low- Wind models. Because the kick must nearly exactly counter 
the momentum lost in the supernova ejecta for our strong supernova models, the systems 
that are more likely to remain bound for these NS binaries are those where the countering 
effects lead to small proper motions. None of the kick velocities are high enough to escape 
their host galaxy's gravitational potential unless the galaxy is extremely small. 

Figure 8b shows the merger times for the same model as Fig. 8a. The peaks below 
100 years arises from binaries following the "Brown Merger" scenario or the "MS Merger" 
scenario (Scenarios II and III)*^. For most of these systems, the stars merge within a few 
months of the collapse of the primary. For Scenario I, the merger rates are determined by 
the remaining main-sequence lifetime of the companion. Note that the low wind (small mass) 
subclass dominates the long merger times, but nearly all of the subclasses have the full range 
of merger times. Figure 8c shows the distribution in distance traveled before outburst for 
the He-merger population. This distribution is determined by multiplying the velocity of the 
binary times its merger time. Here again, the <lpc bin is made up entirely of the "Brown" 
scenario and the main-sequence kick scenario. 

In Figures 8d-f, we show the same set of data for the velocity, merger time, and distance 
distributions as a function of fractions relative to the total number in a subclass instead of 
total number. Because the rate of these mergers is still very difficult to determine due to a 
number of uncertainties (e.g. binary mass fraction and initial mass function), such relative 
abundances are most appropriate for comparison to observations. From these plots we can 
see that in the total He-merger population, less than 20% occur within 1 pc of their initial 
formation site. This is essentially a statement that less than 20% of He-merger systems are 
formed in the Brown and MS Merger scenarios. The moment of inertia subclass produces 
the most extended distribution of mergers. 

To demonstrate the effect of the kick distribution, we compare to a single Maxwellian 
kick distribution. Figure 9a-c shows the fractional distributions of velocity, merger time, 
and distance travelled for our Max300 mode. The scales have been kept identical to those 
used in figures 8d-f and these simulations can be compared directly. Note the lack of the 



"^Thc peak between 10-100 y is produced by the Brown scenario in systems where the hehum cores do not 
go through a common envelope phase but are kicked in such a manner to merger - very rare indeed 
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double peak in the velocity distribution. Figures 10 and 11 show the velocity, merger time, 
and distance travelled distributions for a series of models, but focusing only on the total 
He- merger population. Figure 10 focuses on different mass transfer and binary parameters - 
the two most important being stellar radii and stellar mass loss. Figure 11 shows the results 
of a series of models using single Maxwellian kick distributions of varying strengths. 

Tables 5,6 and 7 summarize the results of all of the simulations. Table 5 shows the 
merger rates assuming that a massive binary is formed every lOOy. In general, the standard 
He-merger scenario dominates > 75% of all the He-mergers formed. We confirm the results of 
Fryer et al. (1999) that the main sequence merger rate is very low (< 5 — 10% of all mergers). 
Given the high uncertainty in the rate, we can not rule out any scenario or subclass, or claim 
any subclass as the sole formation scenario or subclass based on the GRB rate, but both 
the main sequence merger scenario and the high moment of inertia subclass definitely have 
rates on the low end to be the sole path for GRB formation. Table 6 lists velocity and 
merger time information. High velocity systems are rare, so any observations that argue 
for a high velocity in the merging system would place strong constraints on some of the 
population synthesis uncertainties. The strongest constraint at the moment may well be the 
number of systems that merge within a wind background versus systems in a more constant 
medium. The effective observational constraint on this is the distance these systems travel 
before merging (Table 7). Systems with stars that initially collapse to form black holes tend 
to move faster but merge quicker, and are likely to be closer to their formation sites. As the 
observational implications get stronger, it is possible that we will be able to home in on the 
true GRB progenitor subclass. 

Some final summarizing comments on the nature of these populations: 

• Roughly 20% (15% if we exclude main-sequence mergers as they may not be dense 
enough to produce GRBs) of all He-mergers will occur with the wind-driven shroud 
produced by the collapse of the star forming the compact remnant for the GRB. These 
Bursts will not only occur within the wind, but also within a young supernova remnant 
(< lOy old). 



• 



The rest of the systems produce bursts ~l-50 million years after the collapse of their 
massive component, roughly 10-100 million years after the initial burst of star forma- 
tion. 

Although some bursts can occur 1 kpc away from their initial star formation, the bulk 
occur within 100 pc of their formation sites. 

What do these results tell us about the environments around He-Merger GRBs? Although 



• 
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the star forming the compact object in these systems may have had a Wolf-Rayet wind profile, 
because of the delay in the explosion, the environment surrounding the system during the 
GRB outburst need not have this r~^ density profile. If all of the progenitors we discussed 
here make GRBs, 80% of He-mergers will not be characterized by a strong Wolf-Rayet wind 
environment. If the companion to the compact remnant exhibits strong winds, its wind profile 
will dominate the surrounding medium. But, in our simulations, 60-90% of the companions 
have masses below 150 and hence will not have winds that strongly affect the surrounding 
medium. In comparing to GRB afterglows, the surrounding media of these systems (80% 
X 60-90% = 50-70% of all systems) will exhibit constant density profiles. But take these 
fractions with some caution - we do not know what subset of all He-mergers will actually 
produce GRBs, and depending upon this subset, the relative rates will change. 



4. Conclusion 

We have studied the environments around coUapsar-like long-duration gamma-ray bursts. 
We have focused on two classes of progenitors that have been proposed to explain these GRB 
explosions with their associated supernovae: the collapse of a massive star and the merger 
of a compact object with its stellar companion. The progenitor class leading to massive star 
collapse is characterized by progenitors whose mass loss prior to collapse is extremely high. 
The environment surrounding these progenitors is defined by this mass loss rate. The general 
assumption (e.g. Chevalier et al. 2004) is that the density profile follows an r~^ profile with 
a density proportional to a constant mass loss. But the mass loss for these massive stars is 
anything but constant. This variability leads to modifications to the r^^ profile. In addition, 
the mass loss can be significantly asymmetric, leading to additional wrinkles in the profile of 
the surrounding environment. In general, these modifications are not sufficient to resolve the 
discrepancies between predictions of this progenitor scenario and the observations. Recent 
work studying mass loss in GRB progenitors has indicated that these environment observa- 
tions contain a wealth of information about the bursts (van Marie et al. 2005, Eldridge et 
al. 2006). 

Although further study (both theoretical and observational) is required, we are now 
pushed to some extreme predictions in order for the classic massive star progenitor scenario 
for GRBs to produce roughly half of all GRB explosions within constant density environ- 
ments. In addition, the low densities predicted by some GRB progenitors will be much 
harder to produce with massive star progenitors. This, in effect, is strike two for the mas- 
sive star collapse progenitor for long-duration GRBs - the first strike being the difficulty in 
getting enough angular momentum to form black hole accretion disks in these progenitors 
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(e.g. Petrovic et al. 2005, Fryer & Heger 2005). A possible solution may still exist if GRB 
progenitors are low-mass (20-25 Mq), low-metallicity stars. The resultant Wolf-Rayet winds 
for such stars could be less than IO^'^MqI/"^. 

The other class of progenitor, the He-merger model, invokes the merger of a compact 
remnant from a massive star with its companion. Although this progenitor model also 
begins with a massive star with strong winds, the system can travel up to Ikpc from this 
wind environment prior to the GRB explosion. We described 3 different formation scenarios: 
the classic "He-merger" model produced when the companion envelops its compact remnant 
companion, the Brown merger scenario requiring nearly equal-massed binary stars, and the 
main-sequence kick scenario. The relative rates of these three scenarios are roughly 80:15:5 
respectively. Essentially all of the classic He-merger systems explode more than 10 pc from 
their initial formation site. For these models, the surrounding environments are not set 
by any wind profile (unless the companion star also has a strong wind). In contrast, the 
Brown and main-sequence merger scenarios occur within the wind profile of the compact 
remnant forming star. Indeed, the GRB from these scenarios will occur within 10-100 y 
of the collapse of the primary (and its possible explosion). Assuming that these scenarios 
produce all long-duration GRBs, we predict at least 20% of all GRBs should have wind 
profiles. The remaining 80% will exhibit weaker wind environments (not from Wolf-Rayet 
stars) of the NS/BH companion, 60-90% of which are stars with masses below 15 M0. This 
means that, if all He-mergers produce GRBs, we expect that 50-70% of all GRBs should 
occur in ISM-like surroundings. Roughly 20-30% of these bursts will occur more than 100 pc 
from their initial formation sites. 

It appears that the He-merger class of long-duration progenitors fits the constant density 
surroundings better than the stellar collapse progenitors. This is not surprising. Fryer & 
Woosley (1998) proposed this model, in part, to explain this observation. However, this 
model still has one strike against it: these mergers tend to have too much angular momentum. 
Note also that these progenitors predict significant hydrogen in an excretion disk not too 
far from the collapsing star. Although one would not expect hydrogen lines to be strong in 
the supernova explosion accompanying the GRB from this He-merger progenitor, ultimately 
evidence of the supernova shock interacting this disk may produce some hydrogen emission 
akin to the hydrogen line emission observed in the Type la supernova 1999ee (Mazzali et al. 
2005). 

If He-merers are the only long-duration GRB progenitor, these constraints will probably 
require that only a subset of all the possible He-merger models actually produce GRBs, but 
much more detailed calculations (both of population synthesis and of common envelope 
mergers) are required to constrain this subset. Observations combined with population 
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synthesis studies may provide clues to the true subclass of He-mergers responsible for the 
formation of GRBs. 

If observations can constrain the subclass of He-mergers that produce GRBs, we may 
also be able to constrain some of the uncertainties currently plaguing population synthesis 
calculations. We can use these comparisons to help validate our population synthesis codes 
for use in other binary calculations: e.g. the prediction of populations of X-ray binaries. 
These calculations also ultimately teach us about stellar evolution. The most clear-cut 
examples of uncertain parameters in this study are stellar radii and stellar winds. In this 
way, studying the environments around GRBs not only teach us about GRBs, but also the 
objects that make GRBs. 

Alternatively, both progenitor scenarios may produce GRBS: the He-merger model 
would be responsible for producing all p =constant ambient medium bursts along with 
some of the p oc r~^ bursts; the classic massive star progenitor will account for a fraction of 
the p oc r~^ wind-blown bubble bursts. 
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Table 1. Wind Parameter Study 



Simulation 


juOy (]\T0.25Myr\ 
zone \ zone / 




M 


^ISM 


''^wind 


^termination 


''"outer 


Name 


(10*^ zones) 


{M^y-') 


(cm-3) 


(lOOOkms-1) 


V termination j 


V outer j 


plO 


6.2(8.7) 




10-5 


10 




2.4 


12.0 


plOO 


7.5(10.0) 




10-5 


102 




1.0 


7.7 


plOOcooP 


7.5(10.0) 




10-5 


102 




1.0 


7.7 


plOOO 


9.3(11.8) 




10-5 


103 




0.42 


5.1 


plOOOO 


12.1(14.6) 




10-5 


1.8 X 10^ 




0.12 


3.2 


Mx 10 


6.2(8.7) 




10-^ 


10 




4.8 


18.4 


M/10 


6.2(8.7) 




10-6 


10 




1.1 


7.4 


'i'wind X 2 


6.2(8.7) 




10-5 


10 


2 


2.0 


16.0 


t'wind/S 


6.2(8.7) 




10-5 


10 


2 


2.0 


13.1 


16 M© 


2.3(2.4) 


10- 


-7 - 10-5 


100 


2-3 


r^4:5^ 


N/A 


23 M© 


2.3(5.0) 


10- 


-7 - 10-1 


100 


0.2-3.5 


-2^ 


N/A 


40 Mq 


2.3(2.8) 


10- 


-6 _ 10-1 


100 


0.3-4.0 


-2^ 


N/A 



'^This simulation uses the cooling routine outlined in Fryer et al. (2006) based on the cooling 
rates of Sutherland & Dopita (1993). 

^The termination shock is less well defined for these models. Instead we give the radius at which 
a roughly p oc r^^ density profile ends. There is certainly structure within this radius, but probably 
not a sufficient deviation to change the conclusions of the afterglow constraints. 
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Table 2. Merger Scenarios 



Name 



Merger Mechanism'^ Companion'^ Merger Time'^ 



Distance'^ 



He-Merger*" 


CE 


Evolved 


10^-lO^y 


~ 10 lOOOpc 


Brown Merger^ 


Kick 


Evolved 


< 100 y 


<0.1pc 


MS Merger'' 


Kick 


Main-Sequence 


< 100 y 


<0.1pc 



'^The NS/BH merges with its companion either by being kicked into it during the super- 
nova explosion (Kick) or by merging when the secondary expands off the main-sequence 
surrounding the NS/BH in a common envelope (CE). 

^We differentiate the companion based on whether it has an evolved helium core or the 
unevolved core of a main-sequence star. The density in the evolved helium core can be much 
higher than that of a main-sequence star. 

'"Time between the collapse of the primary (forming the compact remnant) and the merger 
of this compact remnant with the companion star. 

■^Distance traveled by the system after the primary collapse. 

^Fryer & Woosley (1998) 

^Fryer et al. (1999) 



Table 3. Subclasses 



Name 



Constraint 



Trend in Progenitors 



NS Collapse Afprimary > 2OM0 Low-Mass Primary 

BH Collapse ^Primary < 2OM0 High-Mass Primary 

Low- Wind ^secondary < ISM© Low-Mass Secondary 

Angular Momentum JneStar > /compactBinary High-Mass Secondary 
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Table 4. Simulation Properties 



Simulation 


aiMp'' 


^^kick"" 


o:mt,Pmt'^ 


«ce'^ 




Other 


Standard 


2.35 


FBB 


1.0,0.8 


0.5 






StanCEO.5 


2.35 


FBB 


0.5,0.5 


0.5 






StanCEl.O 


2.35 


FBB 


0.5,0.5 


1.0 






StanCEO.2 


2.35 


FBB 


0.5,0.5 


0.2 






IMF2.7 


2.7 


FBB 


1.0,0.8 


0.5 






IMF2.7CE0.5 


2.7 


FBB 


0.5,0.5 


0.5 






IMF2.7CE1.0 


2.7 


FBB 


0.5,0.5 


1.0 






IMF2.7CE0.2 


2.7 


FBB 


0.5,0.5 


0.2 






IMF2.7MT1.0 


2.7 


FBB 


NA,1.0 


0.5 






FBBHighWind 


2.7 


FBB 


NA,1.0 


0.5 


^^Wind 


= 10 X Mf- 


FBBLowRad 


2.7 


FBB 


NA,1.0 


0.5 


-Rstar 


= 0.1 X Rill- 


FBBHighRad 


2.7 


FBB 


NA,1.0 


0.5 


-Rstar 


= 4 X Rill- 


Max50 


2.7 


50 


0.5,0.5 


0.5 






MaxlOO 


2.7 


100 


0.5,0.5 


0.5 






MaxlSO 


2.7 


150 


0.5,0.5 


0.5 






Max200 


2.7 


200 


0.5,0.5 


0.5 






Max300 


2.7 


300 


0.5,0.5 


0.5 






Max400 


2.7 


400 


0.5,0.5 


0.5 






MaxSOO 


2.7 


500 


0.5,0.5 


0.5 






Max300IMF 


2.35 


300 


0.5,0.5 


0.5 






MaxHighWind 


2.7 


300 


NA,1.0 


0.5 


^Wind 


= 10 X M|^- 



'^For our models, we sample both the primary and the secondary from an 
initial mass function. We set the minimum primary mass to 9 Mq. The 
secondary is allowed to be as low as 1 Mq and we insure that its mass is less 
than the primary. 

'^We use two kick distributions: a bimodal Maxwellian with mean veloc- 
ities at 50 and 500 km s^^ based on the bimodal distribution from Fryer, 
Burrows, & Benz (1998) denoted FBB, and a single Maxwellian denoted by 
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the mean velocity. 

'^Mass transfer parameters: (3ut is the fraction of the expanding star's 
envelope that is accreted onto the companion (the rest is assumed to be 
lost from the system), auT is the fraction of the specific angular momentum 
carried away by the mass that is lost. These parameters both determine the 
mass accreted by the companion and the orbital separation of the binary. 

'^If the expanding star envelops its companion, the companion will spiral 
in through the envelope. acE is a parameter designed to represent how 
efficiently the orbital energy ejects the envelope. It includes assumptions 
both about the binding energy of the expanding envelope and the amount 
of orbital energy injected during the inspiral. If this value is low, the final 
separation after a common envelope phase will also be small. 
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Table 5. Merger Rates (Myr 



Simulation 


Formation Scenario 






Subclass 






Helium 


Brown 


MS 


NS 


BH 


Low- Wind 


Ang. Mom. 


Standard 


150 


26 


1 


130 


60 


150 


7.8 


StanCEO.5 


100 


25 


7.3 


74 


60 


95 


7.8 


StanCEl.O 


99 


23 


7.2 


71 


58 


96 


8.6 


StanCEO.2 


100 


28 


7.2 


77 


59 


91 


5.6 


IMF2.7 


160 


19 


4.5 


140 


40 


160 


5.0 


IMF2.7CE0.5 


103 


19 


4.5 


82 


45 


98 


5.8 


IMF2.7CE1.0 


96 


18 


4.4 


68 


40 


95 


5.6 


IMF2.7CE0.2 


100 


20 


4.3 


84 


40 


93 


3.6 


IMF2.7MT1.0 


68 


19 


4.5 


51 


41 


65 


5.0 


FBBHighWind 


160 


4.3 





149 


12 


149 


0.24 


FBBLowRad 


2.7 


38 





21 


20 


6.6 


1.1 


FBBHighRad 


110 


5.4 


19 


100 


35 


120 


4.4 


Max50 


0.13 


0.026 


0.0064 


0.11 


0.056 


0.13 


0.0066 


MaxlOO 


120 


30 


7.5 


96 


65 


120 


8.1 


Maxl50 


72 


22 


6.7 


52 


49 


69 


6.3 


Max200 


46 


16 


5.4 


32 


36 


44 


4.7 


Max300 


18 


7.9 


3.0 


12 


17 


17 


2.4 


Max400 


5.7 


3.5 


1.3 


3.4 


7.1 


5.1 


1.1 


MaxSOO 


1.5 


1.4 


0.46 


0.78 


2.5 


1.3 


0.38 


MaxSOOIMF 


21 


11 


4.9 


11 


26 


29 


3.8 


MaxHighWind 


19 


1.4 


0.0 


19 


1.3 


18 


0.04 
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Table 6. Velocity and Merger Times 



Simulation 


Velocity > 50km s 


-^ (> 100km f 


.-l)a 


Mean Merger Time^ 




NS 


BH 


Low-Wind 


Ang. Mom. 


(1 — 0" deviation) 


Standard 


0.31(0.04) 


5.9(1.4) 


5.0(0.47) 


0.549(0.019) 


9.9(5.6) 


StanCEO.5 


0.54(0.062) 


4.1(0.53) 


4.9(0.46) 


0.83(0.030) 


10.4(7.9) 


StanCEl.O 


0.54(0.065) 


2.9(0.15) 


4.4(0.41) 


0.83(0.030) 


10.6(8.1) 


StanCEO.2 


0.48(0.062) 


7.7(1.7) 


6.8(0.68) 


0.82(0.033) 


9.9(7.3) 


IMF2.7 


0.26(0.030) 


5.3(1.1) 


5.0(0.41) 


0.36(0.011) 


10.4(5.5) 


IMF2.7CE0.5 


0.46(0.056) 


4.1(0.42) 


5.4(0.39) 


0.70(0.018) 


11.0(7.9) 


IMF2.7CE1.0 


0.49(0.061) 


2.7(0.12) 


4.3(0.34) 


0.61(0.016) 


11.3(7.9) 


IMF2.7CE0.2 


0.42(0.055) 


7.0(1.5) 


6.6(0.60) 


0.57(0.018) 


10.4(7.1) 


IMF2.7MT1.0 


0.74(0.088) 


3.0(0.13) 


4.8(0.40) 


0.87(0.024) 


10.7(8.8) 


FBBHighWind 


0.077(0.0024) 


0(0) 


0(0) 


0(0) 


10.4(5.4) 


FBBLowRad 


24(14) 


24(10) 


12(2.0) 


3.6(0.63) 


9.5(9.9) 


FBBHighRad 


0.009(0) 


2.4(0.096) 


3.2(0.11) 


0.25(0.0018) 


10.1,4.8 


Max50 


0(0) 


0.72(0) 


3.0(0) 


0.32(0) 


11.1,(7.0) 


MaxlOO 


0.42(0.0023) 


2.3(0.0046) 


2.9(0.0074) 


0.56(0.0034) 


11.0(8.1) 


MaxlSO 


1.1(0.0054) 


4.5(0.015) 


5.6(0.015) 


1.4(0.0015) 


11.0(8.6) 


Max200 


3.1(0.016) 


11(0.036) 


13(0.043) 


3.7(0.0041) 


11.0(8.8) 


Max300 


11(0.10) 


35(0.27) 


40(0.21) 


15(0.024) 


11.0(9.2) 


Max400 


23(1.2) 


49(1.7) 


55(1.1) 


24(0.19) 


10.9(9.8) 


MaxSOO 


44(9.6) 


65(12) 


69(7.8) 


39(1.8) 


11.0(10.7) 


MaxSOOIMF 


13(0.13) 


37(0.28) 


41(0.22) 


19(0.034) 


10.1(9.3) 


MaxHighWind 


2.4(0.0052) 


0.015(0) 


0.54(0) 


0(0) 


10.4(5.4) 



'^Percentage of systems in each subclass that have velocities greater than > 50km s ^ (> 
100km s~^) 

'^Mean merger time in Myr (1 — a deviation in this merger time assuming Gaussian distribution). 
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Table 7. Distances Travelled between collapse and Merger 



Simulation 


Distances > 


> lOpc (> 10( 


)pc)^ 




NS 


BH 


Low- Wind 


Ang. Mom. 


Standard 


88(18) 


53(22) 


76(34) 


91(21) 


StanCEO.5 


84(15) 


53(23) 


75(35) 


90(22) 


StanCEl.O 


83(15) 


58(24) 


76(34) 


90(23) 


StanCEO.2 


86(14) 


44(17) 


71(36) 


89(19) 


IMF2.7 


89(18) 


57(26) 


77(38) 


92(22) 


IMF2.7CE0.5 


85(16) 


60(29) 


79(42) 


91(23) 


IMF2.7CE1.0 


84(16) 


62(29) 


68(39) 


91(24) 


IMF2.7CE0.2 


87(15) 


48(19) 


73(40) 


90(19) 


IMF2.7MT1.0 


82(13) 


55(27) 


76(91) 


91(23) 


FBBHighWind 


90(19) 


75(1.7) 


28(2.4) 


92(19) 


FBBLowRad 


0.54(0.20) 


12(6.6) 


92(52) 


27(15) 


FBBHighRad 


84(6.6) 


33(14) 


35(18) 


79(8.7) 


Max50 


87(17) 


51(21) 


58(24) 


91(22) 


MaxlOO 


84(18) 


58(30) 


79(43) 


92(27) 


MaxlSO 


79(17) 


57(34) 


77(49) 


91(31) 


Max200 


76(15) 


56(37) 


77(52) 


90(34) 


Max300 


69(15) 


54(39) 


75(56) 


89(39) 


Max400 


57(13) 


50(39) 


74(60) 


88(47) 


MaxSOO 


36(8.4) 


45(38) 


76(66) 


85(57) 


MaxSOOIMF 


66(14) 


50(34) 


73(52) 


87(40) 


MaxHighWind 


87(23) 


79(2.3) 


26(3.3) 


93(24) 



^Percentage of systems in each subclass that have distances 
greater than > 10 pc (> 100 pc) 



100 F 
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Fig. 1. — Density vs. radius for constant wind models. The top panel shows models with 
mass loss rates of 10~^ Mq|/~^ and wind velocities of 1,000 km s~^ for a range of densities in 
the ambient interstellar medium. The wind extends well into the interstellar medium (beyond 
10 pc for ambient densities of 10cm~^). For afterglow observations, the most critical feature 
is the size of the free-streaming region. It is in this region that the density follows a p oc r~^ 
dependence on distance. This region is smaller for higher interstellar densities. Sound waves 
move back and forth through the shocked wind/ISM region, causing density perturbations, 
both real and numerical. The wiggles in the shocked region of our highest- density cases are 
numerical in origin. The dotted line in this top panel shows the result from our simulation 
including cooling (see Table 1 for details). The bottom panel studies the dependencies on 
the mass loss rate and the wind velocity around a model assuming a 10~^ MqH'^ mass 
loss rate, wind velocity of 1,000 km s~^ and ambient density of 10cm~^. The solid lines 
compare different mass loss rates: dark (IO'^'MqU"^), light (10~^ Mq-u^^). The dotted lines 
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Time (y) 



Fig. 2. — Mass loss rate (top) and wind velocity (bottom) time histories for three stellar 
models: 16 and 23 Mq binary systems and a 40 Mq single star. Note that the mass loss is 
far from constant with some pulses lasting for less than 100,000 y. The largest mass outflows 
often are the lowest velocity flows and this can lead to shocks as the fast moving material 
following the peak outflows catches these peaks. 
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Fig. 3. — Density versus radius for the wind bubbles produced by our 3 stellar models: 40 Mq 
star (solid line), 23 M© star (dotted), 16 M© star (dashed hue). The effects of the variable 
mass loss can be seen (for instance, see the density peak in the 23 M© star at roughly 3pc), 
but the density profiles still follow a roughly p oc r^^ profile even beyond the true termination 
shock of the free-streaming wind and far enough out that it should dominate what we see 
in the GRB afterglow observations. The points show the density profile along the z-axis in 
our 3-dimensional model. This model was for the 23 Mq binary star assuming zero density 
ambient medium except for a circumstellar disk caused by the binary interaction (see Fig. 
4). 
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Fig. 4. — A 3-D SPH model of wind from a massive star interacting with a wedge of dense 
material ejected during the common envelope phase of a binary system. The location of each 
SPH particle is marked by a colored point, and the color indicates the number density of 
gas at that location. The wind-producing star is located at the center of the simulation box; 
the wedge of material ejected during the common envelope phase is oriented perpendicular 
to the z (vertical) axis. Wind particles are injected at a radius of 0.01 pc, the inner edge of 
the dense wedge is located at 0.3 pc, and particles are removed from the simulation when 
they travel more than 3 pc from the origin. The density of gas in the wedge is much larger 
than in the subsequently-ejected wind, so the wind is obstructed by the wedge and escapes 
primarily along the z axis. 
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Scenario I: Standard (CE evolution) 




Orbital Separation's 1 AU 



Primary Evolves off Main Sequence 
Primary Expands 



( Red Giant |j 



Roche Lobe Overflow 
Mass Transfer, Possible 
Common Envelope 



Collapse - Formation of BH or NS 



({ Red Giant j 



Secondary Evolves 
Off Main Sequence 
X-ray Binary Phase 



Common Envelope Phase 
Orbital Separation Shrinks 
Hypercritical Accretion 



Stars 
Merge 





Compact object Accretes to become BH 

Fig. 5. — Standard formation scenario for He- merger progenitors. Two stars in a close 
binary undergo at least 1 common envelope phase, and often 2 such phases. The first, not 
required, occurs when the primary (more massive star) evolves off the main sequence. After 
the common envelope phase, the primary collapses to form a neutron star or black hole. If 
the system remains bound after this collapse and its orbit is close enough, the binary will 
go through a second common envelope phase when the companion to the compact remnant 
evolves off the main sequence. If the compact object merges with its companion, a Helium- 
merger is formed. 
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Scenario II: "Brown" Scenario 

Initial Conditions: 

Mp> MsN 

Ms^0.95Mp 

Orbital Separation $ 1 AU 

Primary Evolves off Main Sequence 
Primary Expands 





Common Envelope Phase 



Secondary Evolves Off Main Sequence 




Double Helium Star 
Binary 



Envelope Ejected - Primary Collapse 



NS/BH Kicked into 
Helium Companion 




Fig. 6. — Secondary formation scenario of helium-mergers. Two nearly equal mass stars in 
a close binary evolve off the main sequence together - that is, the secondary evolves off the 
main sequence before the primary collapses. The "joint" common envelope phase from this 
evolution produces two helium stars. When one star collapses, its compact remnant may be 
kicked into the helium companion, producing a helium-merger. 
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Scenario III: Kick Progenitor 





Orbital Separation^ 1 AU 



Primary Evolves off Main Sequence 

Primary Expands 

Roche Lobe Overflow 
Common Envelope Evolutioi 



Collapse - Formation of BH or NS 



NS,BH kicked into Main Sequence Companion 
Stars Merge 




Fig. 7. — Alternate scenario for producing merged systems. In this close binary system, the 
collapse of the primary star kicks the compact remnant into its main-sequence companion. 
Although not strictly a Helium-merger because the companion is still on the main sequence 
(and hence not quite as dense and less likely to be an ideal GRB progenitor), we include it 
in our study. 
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Fig. 8. — A series of plots showing the population for all our He- merger formation scenarios 
(either total number or fraction) for our Standard model versus system velocity, merger 
time or distance travelled. The different curves denote different subclasses and have the 
same definitions for each plot: dark solid - all systems, dark dotted - angular momentum, 
light solid - BH Collapse, light dotted - NS Collapse, and dashed - Low- Wind (see Table 
3 for details), a) Number of binaries per Ikms^^ bin. The peak in low velocity systems 
only occurs in the neutron star compact object and small companion mass subclasses. The 
velocities are only high enough to escape the host galaxy gravitational potential if the host is 
an extremely low-mass galaxy. Note that the peak in low-velocity systems is dominated by 
the NS Collapse and Low- Wind models. It is difficult for those systems that form neutron 
stars to remain bound, requiring the kick to counter the momentum lost in the supernova 
ejecta. The systems that are more likely to remain bound for these NS binaries are those 
where the countering effects lead to small proper motions, b) Number of binaries per log 
Time bin (900 bins in plot). The large spike in binaries with merger times < lOy arises from 
the Brown and MS Merger formation scenarios. These scenarios only produce immediate 
mergers. The bulk of long merger times are produced by those mergers with low wind (and 
hence low-mass) companions, c) Number of binaries per 1 pc bin. The spike at low distances 
corresponds to the Brown and MS-merger formation scenarios where the merger time is 
nearly immediate (see figure 8b). Even for the primary He-merger scenario, most mergers 
occur within 100 pc of ther formation site. A small fraction occur beyond Ikpc (Table 6). 
d-f) Same as figures a-c) but showing the relative fractional distribution (Number in a bin 
divided by the total number of a given subclass). 
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Fig. 9. — Relative fractional distributions for model Max300. These 3 figures are identical 
in scale to figures 8d-f and can be compared directly to the results of the standard model in 
figure 8. Note that because the kick in this model is single-peaked, the double peak in the 
velocity distribution does not occur in this model. 
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Fig. 10. — Same as Figure 9 showing just the total He-merger distributions for a range 
of population synthesis progenitors: black solid (Standard), red solid (FBBLowRad), red 
dotted (FBBHighWind), green sohd (StanCEl.O), green dotted (StanCEO.2), blue solid 
(IMF2.7MT1.0), blue dotted (IMF2.7CE0.2), cyan sohd (IMF2.7CE0.5). The largest varia- 
tions are caused by those models where the stellar radii or stellar mass loss was modified. 
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Fig. 11. — Same as Figure 9 showing just the total He-merger distributions as a function of 
kick magnitude (using single peaked maxwellian distributions): blue solid (MaxlOO), cyan 
solid (MAX150), green sohd (MAX200), blue dotted (MAX300), cyan dotted (MAX400), 
green dotted (MAX500), red dotted (MAX300IMF). Note that the kick distribution has a 
much larger effect than changes in the initial mass function. 



-44- 



0.1 






a 




a 
& 


"^v^ 








" 0.01 

1 


; ■■"■.. V 


^^*^S5&^,,:.,, 


.;--■ ..^., 


- 












1 




\\\ \ 






3 0.001 




\ \ H 






1 




\^, 


\ " ■ ''' 




OOOOl 




\; 


1^4 J^ : 

>. In A ;a li. 





Velocity (km s"') 




0.01 



C 

S 0.001 

o 

a 



a 

B 

.2 

o 
Id 



0.0001 



— 1 1 1 1 1 1 1 1 r- 



— 1 1 1 1 1 1 1 1 1 1 1 1 [— 




_l i I L 



J I I L. 



100 



200 300 

Distance (pc) 



400 



500 



Fie. 11. 



